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Very-high energy observations of blazars can be used to constrain the strength of the intergalactic 
magnetic field. A simplifying assumption which is often made is that of a magnetic field of constant 
strength composed by randomly oriented and identical cells. In this paper, we demonstrate that a 
more realistic description of the structure of the intergalactic magnetic field is indeed needed. If such 
a description is adopted, the observational bounds on the field strength are significantly affected in 
the limit of short field correlation lengths: in particular, they acquire a dependence on the magnetic 
field power spectrum. In the case of intergalactic magnetic fields which are generated causally, for 
which the magnetic field large scale spectral index is ns > 2 and even, the observational lower 
bound becomes more constraining by about a factor 3. If instead —3 < ng < —2, the lower bound 
is significantly relaxed. Such magnetic fields with very red spectra can in principle be produced 
during inflation, but remain up to now speculative. 


I. INTRODUCTION 

Very High Energy (VHE) gamma-rays (> 100 GeV) 
from distant extragalactic sources are absorbed by pair 
production on the intergalactic background radiation 
field [I]. Electron-positron pairs produced in this way 
would then inverse Compton scatter off soft ambient pho¬ 
tons, initiating an electromagnetic cascade that might 
be detected in the multi-GeV to TeV energy domain [2|. 
The presence of an InterGalactic Magnetic Field (IGMF) 
affects the development of the electromagnetic cascade by 
deflecting (or even isotropizing) electron- positron pairs 
[HE]. The deflection would also introduce a delay in the 
arrival time of the photons produced in the cascade, with 
respect to the arrival time of the primary VHE photons 
emitted from the source. For bursting sources, the obser¬ 
vation of such a delay would allow to constrain the value 
of the IGMF 0. 

Blazars were the first extragalactic sources to be de¬ 
tected in VHE gamma-rays |4| , and it was soon realized 
that they are ideal targets to search for the cascade ra¬ 
diation (2j. If the IGMF is strong enough to deflect in 
an appreciable manner the electrons and positrons in the 
cascade, then such a secondary radiation would appear 
in the form of a diffuse gamma-ray emission surround¬ 
ing the central point source of primary photons HE]. 
Moreover, the VHE gamma-ray emission from blazars is 
characterized by strong, powerful flares [e.g.EHUg. Due 
to the time delay induced by the IGMF, the cascade ra¬ 
diation would appear as a delayed echo of such flares HU. 

The IGMF can be modeled as a stochastic magnetic 
field, statistically homogeneous and isotropic. Both the 
angular extension © ext and the time delay fdeiay of the 
cascade radiation from blazars are expected to depend 
on the rms strength of the IGMF, but also on its corre¬ 
lation length Xb mug. The IGMF correlation length 
becomes relevant when it is much smaller than the elec¬ 


tron/positron energy loss length for inverse Compton 
scattering, Ab <§C D e : in this case, the charged particle 
would undergo several deflections before losing its energy 
via emission of inverse Compton photons. Its overall de¬ 
flection angle <5 has been estimated considering that many 
stochastic deflections would cause the electron/positron 
to perform a random walk [12] , leading to 

S ^ if X B « D e , (1) 

Kl 

with R b = E e /(eB ) the Larmor radius, E e the elec¬ 
tron/positron energy, and B the IGMF strength. On 
the other hand, if X B D e , the motion of the charged 
particle on the energy loss length scale D e can be approx¬ 
imated as the motion in a homogeneous IGMF, leading 
to a deflection angle which is simply |T2 j. 

S =* ^ if X B > D e . (2) 

Kl 

In the limit of small deflection and neglecting redshift de¬ 
pendences, the angular extension of the cascade radiation 
is directly proportional to the deflection angle, through 
O ext ~ S/t where r = D/D 1 is the optical depth for the 
primary gamma rays, Z) 7 denotes the mean free path of 
the primary photon, and D the distance to the blazar 
H2H33- The time delay, on the other hand, is given by 
tdeiay — D 7 S 2 /2c [12HI4] . It is clear therefore that both 
0 ext and fdeiay inherit, through <5, a dependence on B 
and Xb- thus, a detection of the cascade radiation would 
allow to measure a combination of these two (largely un¬ 
known!) physical quantities [2l l3l fT2l HR] . Since no cas¬ 
cade radiation has been detected yet [TB| [I?] (but see 
the preprint m ). only constraints on the IGMF can be 
obtained (see (76], or [13] and references therein). 

For very weak IGMF the angular extension of the cas¬ 
cade radiation is smaller than the instrumental angular 
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resolution, and thus the contribution from the cascade 
adds up to the primary and point like gamma-ray signal 
from the blazar. For a given value of the IGMF strength, 
-Spsf(Ab), the extension of the cascade emission would 
be identical to the instrument angular resolution. Note 
that, from Eq. Q and from 0 oxt ~ S/t, -Bpsf explicitly 
depends on X B when A B "C D e , as Bpsf oc y/D e /X B 
[El El. For fields significantly stronger than -Bpsf(Ab), 
the cascade emission would become quite extended and 
its contribution to the point like emission would be corre¬ 
spondingly suppressed. By using this argument, Neronov 
and Vovk [T6j obtained a lower limit on the IGMF 
strength of B > 3 x 10” 16 G from the non detection 
by Fermi of the cascade emission from the TeV blazars 
1ES 0299+00 and 1ES 0347-121. The limit is indepen¬ 
dent on As for Ab > D e , while it gradually strengthens 
as ^/D e X B when X B <C D e , c.f. Fig. 2 of TO], 

Similar results were obtained in U2H2U. Such limits 
have been derived after assuming the blazars to be steady 
sources of VHE gamma-rays. However, as pointed out 
in [14] , since blazars are known to be strongly variable 
VHE gamma-ray sources, the assumption of steady flux 
is tenable only over time scales of the order few years, i.e. 
the time window over which the blazars have been actu¬ 
ally monitored with Cherenkov telescopes and proved to 
be steady. If such assumption is made, the lower lim¬ 
its on the IGMF strength relax significantly, because the 
delay time of the cascade emission can easily exceed the 
observing time window. Thus, a more robust lower limit 
on the IGMF strength is B > 10 -18 ...10 -1 ' G, also scal¬ 
ing as l/y/Xs when X B <C D e [Hi [2211251 . The uncer¬ 
tainty in the above quoted value comes mainly from the 
uncertainty in the extragalactic background light [23] . 

The main aim of this paper is to point out that the 
dependence of the deflection angle on the IGMF correla¬ 
tion scale 6 oc V^b (Eq. ([!])), adopted in the literature 
in the limit X B <C D e (see e.g. [E]h is indeed not always 
valid. This affects the constraints on the IGMF. In fact, 
the square root dependence on As follows directly from 
the assumption that the IGMF is composed of chaoti¬ 
cally oriented cells of size Ag < T e inside which the 
IGMF is correlated. Such a configuration does not sat¬ 
isfy the divergence free condition for the magnetic field. 
The square root dependence on As arises because the 
charged particle performs a random walk while it prop¬ 
agates on the path D e , since it crosses several uncorre¬ 
lated cells. This picture applies only if the magnetic field 
spectrum P B (k) oc k nB at large scales (or small wave 
numbers k ) satisfies n B > —2. On the other hand, for a 
very red IGMF spectrum with —3 < n B < —2, the mag¬ 
netic power is not concentrated on a well defined corre¬ 
lation scale but it continuously grows with scales. Con¬ 
sequently, one cannot assume that the charged particle 
moves following a random walk. We will show that, for 
—3 < n B < —2, the dependence of the deflection angle 
on A_b changes from S oc \TXb to S oc A^ B+3 ^“. This di¬ 
rectly introduces a dependence on the IGMF spectrum in 
the lower limit that is derived from gamma-ray cascades. 


In particular, in the limiting case n B —» —3, the lower 
limit becomes independent on X B . We will also show 
that, when the IGMF spectrum is such that n B > —2, 
even if one goes back to the usual behavior S oc \fXvs 
of Eq. <[T|) , the dependence on the IGMF spectral index 
remains, in the form of an overall multiplicative constant 
that depends on n B . 

At scales smaller than the correlation scale X B , it is 
reasonable to assume that the IGMF spectrum is of the 
Kolmogorov type. The spectral index of the IGMF at 
large scales, on the other hand, can directly be con¬ 
nected with the IGMF origin. A magnetic field must 
be divergence-less, and this property directly implies 
that the large-scale power spectrum of a causally cre¬ 
ated IGMF must be blue: more specifically, n B > 2 and 
even [26] 127. This condition must then be satisfied for 
any IGMF that is generated by astrophysical processes, 
and also cosmological processes operating during a phase 
of standard Friedmann expansion of the universe, when 
the causal horizon is finite: be it today, or during the 
radiation or matter dominated eras (for a review of cos¬ 
mological generation processes, see [H HD- 

On the other hand, being connected with causality, this 
constraint on the spectral index does not apply to IGMF 
that root their origin at inflation. During inflation, in 
fact, the causal horizon diverges [29j and magnetic field 
correlations can be built at larger and lager scales. If the 
IGMF is generated by inflation, it can therefore be char¬ 
acterized by a red spectral index. In this case, it becomes 
particularly important to account for the dependence of 
the lower limit from gamma-ray cascades on the magnetic 
field spectral index. Note that in this work we consider 
only non-helical fields, since our aim is only to show how 
the limits from gamma-ray cascades are modified. For 
the case of IGMF with non-zero helicity, and the possi¬ 
bility to detect them, see Refs. ssaEU- 

The remaining of the paper is organized as follows: 
in section II, we briefly revise the properties of the cas¬ 
cade initiated by a VHE photon in the IGM; in section 
III, we define the properties of the stochastic IGMF and 
evaluate the deflection angle; in section IV we apply the 
previously derived expression of the deflection angle to 
the observational bound on the IGMF strength, and dis¬ 
cuss in particular the case of a IGMF generated during 
inflation by the breaking of conformal invariance through 
the coupling of the electromagnetic field with the back¬ 
ground; in section V we conclude. We use CGS units (but 
in the last section, when discussing the IGMF generated 
by inflation, we set c = 1). 


II. ELECTROMAGNETIC RADIATION 
INITIATED BY A VHE PHOTON 

VHE photons do not travel long distances due to pair- 
production interactions with soft photons in the extra- 
galactic background light |pQ . Gamma-ray photons of en¬ 
ergy E l0 interact mainly with soft photons of wavelength 
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A s ~ 1.5(.E 7o /l TeV)/zm (see e.g. 32]). For a blazar lo¬ 
cated at a distance D , the opacity of the pair production 
process is defined as r = D/D 1 , where D 1 is the mean- 
free-path of VHE photons. In the multi-TeV domain the 
following approximate expression can be used [2]: 

» 40k ( ——'j Mpc (3) 

7 V 20 TeV J ’ 

where for simplicity we neglected the effect of the red- 
shift, and where k is a parameter of order unity that 
accounts for the uncertainties in the present knowledge 
of the extragalactic background light (for a review see 

e.g. IS]). 

The average energy of secondary electrons (hereafter 
we refer to both electrons and positrons as electrons) 
produced in this process is E e = E l0 /2. Electrons of 
such energy inverse Compton scatter mainly off photons 
of the cosmic microwave background radiation, of typical 
energy ccmb = 3fcgTcMB ~ 7 x 10~ 4 eV, where k B is the 
Boltzmann constant and Tcmb the temperature of the 
photon field. The cooling length of such electrons is (e.g. 

my 


D P = 


3 m 2 c 4 


Aa T E e u>cMB 


~ 40 


E r. 


20 TeV 


kpc (4) 


while the energy of the up-scattered photons reads (e.g. 

my 


E -=i{^y<™^° A {w^) 2Tev < 5 » 

where wcmb — 0.25 eV/cm 3 is the energy density of the 
cosmic microwave background radiation, <tt the Thomp¬ 
son cross section, and me 2 the electron rest mass en¬ 
ergy. Thus, the secondary gamma-ray emission has to be 
searched in the sub-TeV energy domain. 


III. STOCHASTIC IGMF AND DEFLECTION 
ANGLE 


motivated in the introduction. The IGMF power spec¬ 
trum Pg(/c) can be modeled with two power laws: 


P B {k) = A ■ 


ar **<*> 

(^) ms iffc>fco 


(9) 


where we have introduced a characteristic wavenumber 
k 0 at which the spectrum changes slope. In order to 
make the comparison with previous analyses, we define 
the magnetic field correlation scale as 

= ( 10 ) 

k o 

Therefore, we assume that this is the physical parame¬ 
ter to be eventually constrained by the data, together 
with the IGMF rms strength ■> J(B 2 ). Consequently, the 
normalisation constant A is better rewritten in terms of 
\fW) : using Eqs. Q, ([8| one finds 

« = ^ = d “ 2p * (t) ’ (11) 
which implies then, with Eq. © 

A = ( 27r ) 2 (n B + 3)(m B -3) (B 2 ) 

2 ng+mg k 3 

For the energy density not to diverge, the spectral indexes 
must satisfy the conditions ng > —3, mg > 3. At small 
scales k > kg, the IGMF has been processed by MHD 
turbulence during its evolution from its generation until 
today: it is therefore natural to assume mg = 11/3, 
corresponding to Kolmorogorov turbulence. On the other 
hand, at large scales k < ko the IGMF has not been 
modified by the turbulent cascade: the spectral index ng 
is directly related to the process that generated the IGMF 
and it is one of the relevant parameters characterizing the 
IGMF (see [T3] and references therein). 

The characteristic scale which is usually assigned to 
the stochastic magnetic field as “correlation” scale is the 
integral scale m 


We model the IGMF as a statistically homogeneous 
and isotropic random field with (_Bj(x)) = 0, spatial cor¬ 
relation function 

(J3i(x)Bj(x + y)) = £ij(|y|) (6) 

and energy density corresponding to the squared rms 
strength: pg = ( B 2 )/(8n) = £g(0)/(87r). In Fourier 
space one has then [TJ 

(-Bi(k)-Bj (q)) = (27r) 3 P y (fc)<5(k - q), (7) 

Vijik) = f d 3 ye‘ k y ^( |y|) = |(<% - kikj)P B (k). (8) 

The projector after the second equality in the second 
line follows from the requirement that div B = 0 (see e.g. 
my Note that we only consider non-helical IGMF, as 


A/ = x—y —- f dk k P B (k). (13) 

ott 2 p B Jo 

From Eq. ([9| it appears that this quantity is only well 
defined for ng > — 2. In this case, A j is the scale at 
which most of the magnetic energy is concentrated, and 
it is directly related to the parameter we defined as the 
correlation scale Ag = 27r/fco: 

A/=Ab 5(ng + 2) f ° rnB> " 2 - (14) 

Therefore, when ng > —2 it is natural to assign to the 
IGMF a correlation scale which is simply Ag ~ A/. 

However, the only condition one has to impose on 
ng is ng > —3 for the IGMF energy density not to 
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diverge. Therefore it is (at least in principle) possible 
to have IGMF with very red large scale power spectra, 
—3 < rig < —2. These cannot have a causal origin, but 
can be produced during inflation nn uni- In this case, 
ko = 2 tt/Xb is not a correlation scale in the same sense 
as A i , but it corresponds to the upper cutoff of the IGMF 
spectrum at generation time, which is directly related to 
the process of generation: in general, for an inflationary 
IGMF, it is given by the inverse Hubble scale during in¬ 
flation k 0 ~ H m f. In the following, we keep the definition 
X B = 2Tr/k 0 to be able to describe both cases (ns > —2 
and —3 < ns < —2) with only one parameter. We now 
proceed to demonstrate that, if one accounts for a generic 
large scale spectral index ns, the limit from the gamma- 
ray cascades on the IGMF strength must be modified 
with respect to those presented so far in the literature. 

In order to derive a limit on the IGMF strength from 
gamma-ray observations we need to compute the deflec¬ 
tion angle S of the electrons produced in pair-production 
interactions of VHE photons mm- In order to evaluate 
the dependence of the deflection angle on the statistical 
properties of the IGMF we follow Ref. [33] (for a more 
refined analysis accounting for the three-dimensional de¬ 
viation of the charged particle, see Ref. uni). The basic 
assumption is that the deflection experienced by the elec¬ 
tron in the IGMF is in any case small. Over distances 
smaller than D e the electron energy E e remains approxi¬ 
mately constant, and therefore its motion in the IGMF is 
described by the equation v = ( ec/E e )[\ x B(r)], where 
r denotes the position of the particle. The particle is 
ultra-relativistic with v = cn, and the variation of the 
velocity is v 1 v. Therefore, one can calculate the de¬ 
flection angle over a time interval At as [34] 

ec 2 f t+At 

Av = (5cnj_ = —- / [rixB(r)]dt. (15) 

Re Jt 

Assuming that the particle propagates in the z || n di¬ 
rection almost unchanged by the presence of the IGMF, 
the rms deflection can be expressed as (dz = cdt ) 

(^ 2 ) = (it) 2 x (16) 

f z z+Az dz i j z+Az dz 2 ([n x B(ri)] • [n x B(r 2 )]) ~ 
(±) [ S H - nihj] f z+Az dzi f z+Az dz 2 (B i (ri)B j (r 2 )). 


Eq. § together with the assumption of small deflection 
11"l - r 2 | ~ | z\ - z 2 \ = C gives 


(8 2 )~2&z 



r Az 


UiUj] 


dUiAQ- 


(17) 


This can be calculated in terms of the IGMF power spec¬ 
trum substituting the inverse Fourier transform of Eq. |8]) 
and performing the integral on £ and on the angle dflj, 
(where f d 3 k = f 0 dkk 2 f 4v dflj), to arrive at 


9 / a \ ^ /*°° 

<<2> = (W u) 1 dkkP »<*> * < 18 > 

sin kD e — kD e cos kD e + (kD e ) 2 Si(kD e ) 

. (kD e ) 2 . 


where we set A z — D e to obtain the total deflec¬ 
tion suffered by an electron over the cooling time, and 
Si(fc-D e ) denotes the sinus integral function [ 35 :. The 
above integral can be performed analytically. With 
Rl = v fW)E e le one arrives at a simple expression for 
the rms deflection angle: 




where is an analytic function (given more 

explicitly in the appendix), plotted in Fig. [l] for differ¬ 
ent values of ns- As expected, in the limit A b D ei 
the dependence on hb disappears, while if As <C D e 
the function n (D e /X b,tib) varies with the IGMF large- 
scale spectral index ng. The asymptotic behavior of 

n ( ns ) is such that 



A,/A _i 

V3, V 



7/6 


/ Ab)" B+2 ( A , n B +3 \ n ( Ab \ 3 

\De) ^7r 2 (27r)"B + i nB (n B +2) 2 (3n B +ll) ) ^ 10n B +20 ^ ^ \D C ) 


2 


if X B » D e , 
if A s <C 7/ e ■ 


( 20 ) 


From the upper row of Eq. (201 one sees that, if the cor- back to the usual expression Eq. ([2]): it is as if the particle 
relation scale A b is very large, the deflection angle goes would propagate in a constant IGMF. The small numer- 
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FIG. 1: The function II «flj defined in Eq. (191 for dif¬ 
ferent values of the IGMF large scale spectral index: from 
top to bottom ns —> —3 (black), tib = —2.5 (blue), ub = —2 
(magenta), ub = 0 (green), ns = 2 (red). For large vs. small 
correlation scale, the function reaches the asymptotic behav¬ 
iors given in Eq. p0|. 


ical factor 2/y/3 is due to the fact that we calculate the 
deflection angle a bit more precisely as an integral over 
the charged particle path, while Eq. ([2]) is approximate. 
On the other hand, when the correlation scale is small 1 , 
one does not recover Eq. 0): an explicit dependence on 
ns appears in the deflection angle. If ub > —2, the 


asymptotic behavior given in the second row of Eq. (20) 
is such that the expression in brackets becomes indepen¬ 
dent on A b'- the constant term dominates. Therefore, 
the deflection angle goes back to the usual expression 
s/(S 2 ) oc V D s Xb/Rl , but it st ill differs from Eq. Q 
by a multiplicative constant: \J{tib + 3)/(10 ng + 20). 
For sufficiently red spectra —3 < tib < —2, instead, 


the first term in brackets of the second row of Eq. (20) 


dominates over the constant. The dependence of the de¬ 
flection angle on A b changes with respect to the one 
given in Eq. ([I]), and becomes explicitly dependent on 


n B ■ \/W) °c A 


Gb+3)/2 
B 


Note that the divergence for 


ns = —2 in the analytic expressions of Eq. (201 is only 
apparent. 

This behavior is also seen by comparing directly 


Eq. (19) with the usual deflection angle given in the lit¬ 
erature (Eqs. 0 and 0 ), based on the random walk 
model and therefore here denoted 5 rw : 



& itX ‘ >>D - (21) 


^RW ( D e 

x/M" U’^/h ifAs «D e . 

This quantity is plotted in Fig. [2] 


1 Note that this is usually the case for inflationary generated 
IGMF, for which the upper cutoff ko — H[ n f corresponds to very 
small scales today, even after evolution by MHD processing. 



FIG. 2: The ratio between the deflection angle taken from 
the literature 5 rw (Eqs. 0 and (|2|) the one evaluated here, 
Eq. (211. From bottom to top, ub —► —3 (black), ns = —2.5 
(blueh ns = —2 (magenta), Ub = 0 (green), ns = 2 (red). 
The two angles are the same for large correlation scale As 
(a part from a small numerical factor), but diverge if the 
correlation scale is small, depending on the value of ns- If 
ns > —2, they only differ by a multiplicative constant that 
depends on ns (see main text). The spikes in the figure are 
due to the fact that <5rw is only defined in the limiting cases 
As O’ D e and As D e , and we join them at As = D e . 


IV. IGMF CONSTRAINTS FROM 
GAMMA-RAY CASCADES 


The lower limit on the IGMF strength that has been 
obtained from gamma-ray observations gets modified cor¬ 
respondingly to the change in the deflection angle de¬ 
rived above. In order to illustrate this effect, we consider 
the lower limit from the observations of delayed emission 
which are derived in m in particular, we reproduce the 
result shown in Fig. 14 of this reference. 

Neglecting redshift dependences, the expression of the 
time delay becomes m- 


die lay 


Dry 

2c 


Dry 

2c 


( 5 2 ) 

Pjl\ —tL n 2 

Rl) De 



( 22 ) 


where in the second line we have substituted Eq. (191. 
In the case A b ^ D e , when <5 rw = D £ /Rl, we can 
directly compare the above expression to the one given 
in Eq. (151) of [13], to find 


Idelay 


0.3 


Ery 

0.1 TeV 


(B 2 ) 

(10- 17 G) 2 


A B 
De 


n 2 



yr 


if Xb » De , (23) 


which, choosing = 0.1 TeV and fdeiay = 0.3 yr, di¬ 
rectly translates into the limit 



~ 10” 17 G 




n B 


(24) 
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FIG. 3: Lower bound on the IGMF strength as a function 
of Xb/D e . The grey-shaded region below the dashed curve 
corresponds to the excluded region of Fig. 14 of [X3]. The 
other curves are the lower bounds on y/{B 2 ) for, from bot¬ 
tom to top, ub —► —3 (black), ns = —2.5 (blue), hb = —2 
(magenta), hb = 0 (green), ub = 2 (red). 


This limit is shown in Fig. [3] together with the one from 
Fig. 14 of [T3| , for several values of ns- For large cor¬ 
relation scale, our limit is comparable to the one given 
in m (a part from the negligible small numerical factor 
given in the upper line of Eq. (20)). For small correlation 


scale instead, our limit depends onn^. In particular, if 
ns > — 2 the dependence on As is the same as the usual 
one i J(B 2 ) lim oc \JD e /\g, but our limit is more con¬ 
straining than the one of m by a multiplicative constant 
that depends on ns- This is due to the particular shape 
of the function II: from Eq. (20), one infers the constant 




10 


1-17 



10 tib + 20 


ub + 3 




When —3 < ns < —2, instead, our limit is relaxed with 
respect to the one of Ref. m- In particular, it becomes 
independent on As for a scale invariant IGMF (c.f. again 
Eq. @). 

As mentioned in the introduction, IGMF generated by 
a causal process (anything a part from inflation) have 
blue spectra such that Ub > 2 and even [m w\- In 
this case, we have seen that the usual lower bound from 
gamma-ray cascades should get modified only by a mul¬ 
tiplicative constant depending on ng, but it retains the 


usual dependence on As as yj( B 2 ) hm oc VdT/x B- 
IGMF with spectra ns < 2 could only be generated by 
inflation. For a review of inflationary generation mecha¬ 
nisms, see Refs. Hang. In order to illustrate the effect 
of the new IGMF lower bound found above, we repro¬ 
duce here Fig. 16 of 53, which shows the observation- 
ally testable region for inflationary IGMF in scenarios 
where the IGMF generation is based on the breaking of 
conformal invariance through the coupling of the electro¬ 
magnetic field with the background (via direct coupling 
to the inflaton, or to another scalar held). We reproduce 


the figure for two values of the spectral index: -3 

shown in Fig. [4j and ns = —2 shown in Fig. [5] 

The outcome of the inhationary generation process, i.e. 
the IGMF amplitude, correlation scale, power spectrum, 
depend on the details of generation mechanism. However, 
for scenarios based on the coupling of the electromagnetic 
held to the background, one has in general that the rms 
IGMF amplitude and its correlation scale, redshifted to 
today, are given by 





inf 


( — 

V^reh 


-Hint \ T 0 J \ g 0 J 


(-) 

Vffreh/ 

1/3 


2/3 


0^2 / F re h A / f/ r eh A 


(26) 

(27) 


Here we have inserted the multiplicative factors oq and 
«2 whose amplitude depends on the particular gener¬ 
ation process. Moreover, H ln f is the Hubble scale at 
inflation, To and T re h are the temperatures today and 
at reheating time, and correspondingly go = 3.9 and 
g re h = 106.75 are the number of relativistic degrees of 
freedom. The Hubble scale is such that 3H? nl Mp ~ 
p re h = OsbffrehF r g h , with Mp the reduced Planck mass 
and a s b the Stefan-Boltzmann constant. The initial con¬ 
ditions given in Eqs. (26) and (271 are shown by the green 
line in Figs. [4] and [5] from the maximal reheating temper¬ 
ature of T re h — 10 4 ” GeV down to the minimal one, which 
we set at T re h — 200 GeV as in [13]. To plot y/ (H 2 ) reh 
as a function of A^ h , we have set oq and 0:2 to one in 
Eqs. (26) and (27): this is why we stress that these are 
only “possible” initial conditions. 

The maximal and minimal reheating temperature 
also correspond to a minimal and maximal correlation 
scale, following Eq. (27): setting 0:2 = 1, one has 
A^ h (10 16 GeV) ~ 6 x 1CT 24 Mpc and A# 11 (200 GeV) ~ 
3 x HT 10 Mpc. These are represented by the left, re¬ 
spectively right red-dashed vertical lines in Figs. [4] and 
0 


After generation, the IGMF evolves during the expan¬ 
sion of the universe according to free turbulent decay. For 
a non-helical magnetic field, the evolution is such that the 
final values of the IGMF and its correlation scale evolve 
through selective decay, giving [13] f36] 


\J(B 2 )\ b 2 = const = \J( B2 ) in (Ab) B + , (28) 


with \J{B 2 ) in and Ag the initial conditions. Note that, if 
ub = —2, the selective decay evolution is the same as the 
one of an helical magnetic field (inverse cascade due to 
conservation of helicity). Examples of the selective decay 
evolution are shown by the black arrows in Figs. [4] and [5] 
The endpoints of the evolution by free turbulent decay 
are shown by the blue solid line, and must be such as to 
satisfy the relation [HE?] 

V^> = 10_s iS^ G ' (29) 
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The observationally testable region for inflationary 
IGMF is given by the set of all possible initial conditions, 
i.e. couples of values (yj (7? 2 ) ; „, \g) which, following the 
evolution given in Eq. (28), can fall on the end-line ( f29j ). 
With respect to the result given in Fig. 16 of [T3], here 
we have that in the case ng -)■ —3 the testable region 
opens up to include also the “possible” initial conditions 
given in Eqs. (26) and (27): c.f. the blue-shaded area in 
Fig. @ For ub = —2, on the other hand, these initial 
conditions remain excluded: they are not contained in 
the blue-shaded area of Fig. [5] For intermediate values 
of —3 < ub < —2, the testable region would change as 
shown in Fig. [3] For ng > —2, the testable region dif- 

m 


fers from the one shown in Fig. 16 of 


effect of the overall constant (c.f Eq. (25)): the lower 


only by the 


limit on \J ( B 2 ) becomes a bit more restrictive, and the 
observationally testable region shrinks accordingly. Con¬ 
sequently, testable IGMF with ns > —2 can be produced 
only if the values of aq and a -2 are significantly large than 
one. This can actually happen quite easily: see e.g. the 
model studied in Ref. [SB] and references therein. An¬ 
other possibility is, if the magnetic field gets strongly 
amplified by dynamo action during reheating. The same 
conclusion had been drawn in Ref. [13] . There are cur¬ 
rently no studied mechanisms that can lead to a large 
amplification during reheating, but this is not excluded. 
For example, recently Ref. [39] has studied in details the 
post-inflationary evolution of helical magnetic fields gen¬ 
erated by a pseudoscalar coupling to the inflaton in a 
m 2 (f> 2 inflationary model. 

To summarize, inflationary scenarios that would pro¬ 
vide IGMF with spectral indexes ng > —2 are more con¬ 
strained than previously thought, since the observation- 
ally testable region becomes more restrictive than what 
derived in m ■ On the other hand, scenarios that would 
provide —3 < ub < —2 are less constrained. Note, how¬ 
ever, that valid IGMF generation mechanisms based on 
the breaking of conformal invariance through the cou¬ 
pling of the electromagnetic held with the inflaton al¬ 
ways tend to give IGMF with blue spectra with ns > 1 
m ■ The most widely studied generation mechanism 
which could provide ng -> —3, first proposed in m, 
has a strong-coupling problem [41] . Still, we find it use¬ 
ful to present how the constraints from gamma-ray cas¬ 
cades would change, in case in the future some new, well- 
behaved IGMF generation mechanism is proposed, that 
can lead to red spectral indexes —3 < ns < —2. 


V. CONCLUSIONS 

The non-observation of cascade radiation from blazars 
has been used to put constraints on the strength of the 
IGMF. The IGMF can be described as a stochastic held, 
fully characterized by its rms strength, correlation scale, 
and power-law power spectrum. One expects therefore 
that the constraints inferred from the blazar observa¬ 
tions depend on all these parameters. Indeed, it as been 


1(T 8 
10-ic 
O lO -12 
1(T 14 
10-ie 
10-is 

lO" 24 1(T 19 1(T 14 1(T 9 lO" 4 10 

Ag[Mpc] 


nB —> -3 



FIG. 4: The blue-shaded area corresponds to the observa¬ 
tionally testable region for inflationary generated magnetic 
fields in scenarios with inflaton coupling to the electromag¬ 
netic held. We show here the case tib —► —3. The upper 
and lower gray-shaded regions are excluded respectively by 
the Planck upper bound on the intensity of a primordial mag¬ 
netic held [42] and by the bound given in Eq. (241. The black, 
dashed line is the bound given in Ref. P31- The vertical red 
dashed lines correspond to the minimal and maximal corre¬ 
lation scale, depending on the energy scale of inflation (c.f. 
main text). The green line represents the set of possible ini¬ 
tial conditions given in Eqs. (261 and (27) with aq = a ?2 = 1. 
The black, horizontal arrow represents one example of the 
evolution of the IGMF through freely-decaying MHD. The 
endpoints of the evolution must fall on the blue line given 
by Eq. ( |29| ). We have taken for the energy loss length of the 
electron D e ~ 80 kpc, corresponding to the observed photon 
energy = 0.1 TeV (see section II). 


found that the lower bound on the IGMF rms strength 
depends on the correlation scale in such a way that it 

strengthens as ^ J(B 2 ) lim cx ^/D e /As for small correla¬ 
tion scales \b <C D ei while it is independent of As in the 
opposite regime. However, no dependence on the IGMF 
power spectrum had been identified in previous analy¬ 
ses. The inverse square root dependence on As comes 
directly from the assumption that the IGMF is composed 
of chaotically oriented cells of size A b D e> inside which 
the IGMF is constant and therefore correlated. Here we 
have used a more refined model to evaluate the deflection 
angle experienced by the electron traveling through the 
IGMF, a model that accounts for the properties of the 
IGMF power spectrum. We have found that in general 
the lower limit inferred by blazar observations depends 
also explicitly on the IGMF power spectrum: in partic¬ 
ular, on its large scale spectral index iib ■ If > —2, 
the above mentioned cell-model does effectively describe 
the real situation: in this case, we recover in fact the 

behavior (B 2 ) lim cx \jD e /\ B when \ B -C D e . How¬ 
ever, we have demonstrated that an overall multiplicative 
constant depending on ns must be included in the am- 
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FIG. 5: The same as Fig. |4]but for n B = —2. 


as an example = 10~ 17 G: this value reproduces the 
constraint given in paragraph 5.6.3 and Fig. 16 of [13| . 

For the most relevant scenarios of IGMF generation, 

the usual dependence oc \]D e /X B continues 

to hold, and the lower bounds from blazar observations 
depend mildly on n B : only in the amplitude and not 
in the As-slope. Causally generated IGMF, in fact, are 
all characterized by n B > 2 and even. Very negative 
spectral indexes —3 < n B < —2 can only arise if the 
IGMF has been produced during inflation, and no com¬ 
pletely consistent inflationary model has been identified 
yet which does lead to such red spectra. However, in the 
future some new, well-behaved inflationary generation 
mechanism may be proposed, that can lead to IGMF 
with red spectra. In this case, the present result implies 
that the initial conditions it would have to provide 
to fulfill the lower bounds from gamma-ray cascades 
are significantly less constrained than previously thought. 


plitude, which strengthens the constraint somewhat (c.f. 
Eq. (25l). If instead —3 < n B < —2, the dependence 
on the correlation scale when X B <C D e is changed to 

(H 2 ) 1 im oc (D e /\ B ) ( ' nB+3 ' > / 2 . This means that, if the 


IGMF power spectrum is sufficiently red, one cannot use 
the cell-model since X B is not the typical scale at which 
most the IGMF energy density is concentrated. 

In general, we can conclude that the constraints in¬ 
ferred from the blazar observations must be modified as 
given in Eq. p4| : i.e., denoting Hrw the upper limit from 
previous analyses (in the regime in which it is constant, 
X B D e ), one has 




(30) 


Appendix 

The function H ( n B ) is given by 



with 


where the function n (given in the appendix) 

is plotted in Fig. [I] This equation represents the main 
result of this paper. Note that in Eq. (24) we have taken 


f(x) = sin(x) — £ cos(:r) + x 2 Si(x). 
This can be integrated analytically. 
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